Massive skyrmions in quantum Hall ferromagnets 
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We apply the theory of elasticity to study the effects of skyrmion mass on lattice dynamics in 
quantum Hall systems. We find that massive Skyrme lattices behave like a Wigner crystal in the 
presence of a uniform perpendicular magnetic field. We make a comparison with the microscopic 
Hartree-Fock results to characterize the mass of quantum Hall skyrmions at v = 1 and investigate 
how the low temperature phase of Skyrme lattices may be affected by the skyrmion mass. 

PACS number: 73.40. Hm, 73. 20.Dx 



I. INTRODUCTION 

In the last few years topological spin textures in the 
quantum Hall effect (QHE) have received considerable 
attention |p]-|l^]. The existence of skyrmions can be 
anticipated within the frame work of the Chern-Simon- 
Landau-Ginsburg mean field theory i.e., integrating 
out the charge fluctuations of the composite bosons yield 
an effective model for the Chern-Simon's gauge field. The 
transport properties of the skyrmions can be extracted 
through considering the fluctuating Chern-Simon's gauge 
field, which can be derived by expanding the effective ac- 
tion about its minimum energy solution. One of the lead- 
ing terms of this expansion is the Maxwell action of the 
Chern-Simon's fluctuating field, e.g., (w*/p)^ fcw |J S | 2 . 
to* is the electron effective mass in a host semiconduc- 
tor, and p = 1/(27t£q) is the average electron density at 
v = 1, and £q is the magnetic length. 

Although a skyrmion mass is physically reasonable, in 
the usual minimal field theories |^|,^| skyrmions are con- 
sidered as massless objects. ^From the microscopic point 
of view, e.g., a microscopic Hartree-Fock approximation, 
H[l2| it is not certain if the skyrmions are massive. ]15f| 
The resolution of this ambiguity between the microscopic 
Hartree-Fock approximation and the Landau-Ginsburg- 
Chern-Simon's theory (considered in this paper), is an 
open question. 

In this paper we apply the theory of elasticity [Jl6| to 
investigate the effect of skyrmion mass on the thermo- 
dynamic properties of a skyrmion crystal. We derive the 
collective mode dispersion relations for the Skyrme lat- 
tices at long wave lengths. We then make a comparison 
with the massless microscopic Hartree-Fock calculations 
to reconcile the prediction of the Chern-Simon theory 
with the massless models In addition, we suggest 

how the mass of skyrmions may be characterized within 
the microscopic picture. We also study the stability of 
the Skyrme lattices at low temperature and show that 
the low temperature phases of these Skyrme lattices are 
not affected by including a mass term in the effective ac- 



tion, unless if the mass is sufficiently large. We show the 
mass of skyrmions is suppressed by decreasing the Zee- 
man energy, and indicating that skyrmions are massless 
at zero Zeeman energy. 



II. SKYRMION MASS 

If skyrmions have mass, this will affect a host of prop- 
erties, from their tunneling through a constriction to the 
thermodynamics of a lattice of skyrmions. The rationale 
for their having mass is that in the starting point for 
many calculations, the Chern-Simons Lagrangian of Lee 
and Kane, the electrons have mass. After standard ma- 
nipulations (introduction of a CP 1 field, changing vari- 
ables to m, the local spin field) one has a continuum the- 
ory in which gradients in the spin texture become asso- 
ciated with the charge distribution. If a skyrmion moves 
slowly across the system, this corresponds to the motion 
of one quasiparticle. It seems reasonable that the motion 
of such a texture will involve inertial terms. Questions 
similar in spirit arisen in determining the mass of vortices 

us 

Recalling the duality relation between the topological 
3-current of skyrmions and the Chern-Simon's gauge field 
(TLB, — [y ' l^)(.yi V \d v A\, one might anticipate a sec- 
ond derivative term in time will appear in the effective 
action. After integrating out the statistical gauge field 
one obtains in the limit of low frequency and long wave- 
length: § 



S £ [m] = i^(V(fc)|J s| 2 +y|J 
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where the last three terms are the nonlinear sigma action, 
the Zeeman contribution, and a term in the action that 
guarantees the skyrmions obey fermi statistics. The first 



1 



term is electrostatic, and involves the fourier transform of 
an effective interaction potential V(k) including screen- 
ing by fluctuations in the texture, the second contains the 
kinetic energy and the third reflects that the skyrmion 
sees the original boson as a magnetic flux tube. The con- 
stant a is an odd integer and arises from a Chern-Simons 
term that makes the system fermionic. 

The skyrmion current density is related to the spin 
texture via: 

J\ = s—^Xfiu (<9 M m x d v ia) ■ m 

where the indices run over time and two spatial dimen- 
sions. The zeroth component is the topological charge 
density of the texture, which is proportional to the 
quasi-particle number density, p(r) = Jq- If we con- 
sider a single skyrmion texture that moves uniformly, 
J| = J|(x — vt), then it is straightforward to show the 
relation between the skyrmionic current and their charge 
density satisfy the usual charge conservation law, i.e., 
J s = pv. Then the kinetic term simplifies to 

which yields a transport mass for the skyrmions Mq = 
(m*/p) J d 2 rp 2 (r)/(2w) 2 . 

This mass is derived from kinetic considerations, ft is 
possible that the correct "mass" to calculate will depend 
upon exactly what is being measured in a given experi- 
ment. 

III. COLLECTIVE MODES 

The long range order of the skyrme crystal is deter- 
mined by the repulsive Coulomb interaction, and the 
topological, XY interaction of the hedge-hog fields. An 
antifcrromagnetic ordering between the single skyrmions, 
within a square lattice minimizes the topological interac- 
tion. This can be realized after mapping of the topolog- 
ical hedge-hog fields of the charge one skyrmions onto 
a system of classical dipoles. However, a triangular 
lattice is favored by the Coulomb interaction, similar 
to the Wigner crystals. When the Zeeman energy is 
small enough, the skyrmions can pair up into charge 
two skyrmions and lower the total energy of the lat- 
tice. In contrast to single skyrmions where their topolog- 
ical hedge-hog fields are analog to a system of classical 
dipoles, the charge two skyrmions mimic a system of clas- 
sical quadrapoles. This favors triangular lattice ordering 
of charge two skyrmions, i.e., a bi-skyrmion lattice. 

The low-lying collective modes of a Skyrme lattice con- 
sist of phonons and spin waves. |||To|,[l2| The dispersion 
relation of these collective modes can be obtained by 
adding a dynamical term to the effective Hamiltonian. 



This Hamiltonian has been derived in Ref. Q by a first 
principle calculation of a non-linear a model, assuming 
a specific skyrmion is localized and well separated from 
other skyrmions, i.e., they interact weakly. More pre- 
cisely, we assume that we can divide the 2-dimensional 
configuration space into N regions (N is the number of 
skyrmions) such that a given skyrmion lives in the region 
that other skyrmions are close to their vacuum. This 
condition enables us to linearize the interaction poten- 
tial energy among the isolated-skyrmions. This assump- 
tion may break down if the size of skyrmions, A, becomes 
comparable with the distance between them. In this case 
the next to the linear terms in the potential energy may 
be significant, and one should take them into account. 
Roughly speaking, this happens if R < 2A (R is the sep- 
aration between two skyrmions). Here the skyrmion size 
A is defined as the radius at which the spin lies in the 
XY plane. The relevant potential energy functional of 
the excitations is obtained by introducing the lattice of 
equilibrium positions of the skyrmions Ri a is an initial 
lattice ansatz (a is the cartesian component of the posi- 
tion of the ith skyrmion), the displacement field u Q (R;), 
and the orientation field of the skyrmions 0(R$). The 
result is then 

E[u, 0} = Y, VIR, + u(Ri) - Rj - u(R,-)|) 

+ J (i R < + u ( r o - r j - u ( R j)i) 

<ij> 

x cos(x(R,) + 0(Ri) - x(Rj) - 0(R;)). (2) 

Here the topological hedge-hog interaction for a single- 
skyrmion lattice is given by J(R) = c 2 ^/ {^tt 2 £q)K {kR) 
where g — ge 2 /(2e£o) is the Zeeman energy, g is the 
effective gyromagnetic ratio, k 2 — gj \2ir£f } p s ) 1 p s is 
the spin stiffness. Kq(x) is the modified Bessel func- 
tion. In addition, c is a constant that can be ob- 
tained from the asymptotic form of an isolated skyrmions 
and equals c = 30.4£o- |§ For a bi-skyrmion lattice 
J(R) = -c 2 g 2 /{8ir 3 £%p s )K ( K R) and c = 79^. One 
should note c, and therefore the topological interaction 
between skyrmions J(R), depend on the local form of 
skyrmions. 

A. Phonons 

We find the spectrum of the phonons, using the stan- 
dard technique of expanding the energy functional about 
its minima, i.e., u = 0, 8 = 0, assuming the orientational 
field of skyrmions is frozen out. For the single-skyrmion 
case J(R) is positive, hence Xi ~ Xj — ^ is the lowest 
energy state of the topological XY interaction. For the 
bi-skyrmion case J(R) is negative and Xi ~ Xj = is the 
lowest energy state. An estimate on the total energy of 
skyrmions show the bi-skyrmion configuration is likely if 
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the Zeeman energy is very small (< 10 -5 ). We there- 
fore do not consider this configuration for the rest of this 
paper. Expanding the potential energy in terms of the 
displacement fields (up to quadratic order terms) gives 



E[u] = E, 



classic 



^u^k^k^k), (3) 

keBZ a/3 



with ^classic the classical groundstate energy of the 
Skyrme lattice and £) Q| a(k) is the dynamical matrix 

2?r e 2 

D a g(k) = kakp + dcKfi + X)k a kp 

k ea c 

+ f ik 2 5 a/3 + 1 k x k y (l-5 afJ ) + 0(k 4 )]. (4) 

The first term comes from the electrostatic interactions 
of a Wigner crystal, [[l7| the others arise from standard 
2D elasticity theory. The quantity a c is the area of a 
unit cell, while fi, A, and 7 are the conventional Lame 
coefficients for a square lattice. |l6[ ] 

Since the effective interaction consists of two terms, 
i.e., the direct and exchange interactions, the Lame co- 
efficients can be expressed by means of /i = /io + fii, 
A = Ao + Ai and 7 = 70 + 71. Here and 1 are labels 
associated with the direct and the exchange Coulomb 
energy. Expanding the exchange term (second term) in 
Eq.(||) up to the quadratic order in displacement field, 
and for the single-skyrmionic square lattices with only 
nearest neighbor exchange interactions we find 
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(5a) 



(5b) 



(5c) 



We also find the contribution of the direct interaction to 
the elastic constants in units of (e 2 / eloa c )yJ\l — v\/2n 
are /i = -0.146289, A = -0.536199 and 7o = 
— 1.560224 (we take the nearest neighbor separation be- 
tween the skyrmions R = Iq^J2k /\v — 1 as the square 
lattice constant and therefore a c — R 2 ). Note that 7 = 
for a triangular lattices. 
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FIG. 1. The skyrmion mass 

(Ms = 27rMo/m* = J d 2 rp 2 (r)) is plotted as a function of 
the g-factor. The Landau level filling factor is v and m* is 
the effective mass of electron in the semiconductor. 

The time-derivative parts in the action of the displace- 
ment field u, consists of the Wess-Zumino term and the 
kinetic energy of the skyrmions which are first and second 
derivatives in (imaginary) time, respectively. Combining 
these with our previous results gives the following effec- 
tive Euclidean action for the low energy spectrum of the 
displacement field 



1 f m 

Sefflu}^-^ X! / dr [- ie a/3™w c u Q (k,T)9 T u / 3(k,T) - M 5 a pu* a (\<i, T)d 2 Uf3(k,r) + u* (k, r)L) Q/3 (k)u Q (k, t)] , 

^.R 1, r- T3 V J 



aP keBZ 



(G) 



where uj c is the cyclotron frequency of the electron. The 
dependence of the skyrmion mass on the Zeeman en- 
ergy can be obtained by solving the non-linear differ- 
ential equation for the spin texture of a single skyrmion 
H and is shown in Fig. [j]. Note that unlike the to- 
tal number of spins participating in the skyrmion, the 
skyrmion mass decreases with g. This arises because 
the mass is proportional to the square of the skyrmion 
density, and as g decreases the skyrmion becomes more 
spread out. Expanding the time dependence of the com- 
plex field u(k, t) in terms of the bosonic Matsubara fre- 



quencies, uj n = 2nn/P, leads to the fluctuation matrix 

<j /.x _ ( MqujI + D xx (k) muj c uj n + D xy (k) \ , , 
° nW \ -mu) c u> n + D yx {\£) M Lu 2 n + D vy (k) J ' {l) 

The action in Eq. ^ is similar to the action of a Wigner 
crystals in the presence of a magnetic field. Although this 
is as expected, it should be noted that the magnetic field 
interaction is here exactly recovered by the topological 
Wess-Zumino term. 

The contribution of the zero-point energy of the 
phonons to the total energy of the Skyrme lattices, may 
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be obtained by integrating out the quadratic fluctuations x (oj + (k) — ioj n )], 

in u. It turns out that E = Classic + Ej\ n where 



E Ru = lim iV V \n{[3 2 h 2 (u-(k)+iLj n ) (8) 



0— >oo /3 k ^^ z and w± (k) are the phonon frequencies 



2Mg 



w ±( k ) = ^772 m w c + M (£> ra + £>„„) ± J \m?ujl + M (D XX + D yy )Y - AM 2 (D xx D yy - D xy D yx ) ) , (9) 



which have been obtained by substituting the analytical 
continuation iu) n — > w(k) + iS into the fluctuation de- 
terminant. The dispersion relation of the phonons con- 
sists of a gapped mode, as well as a gapless modes. For 
the former, the gap starts at the cyclotron frequency 
uj(k — ► 0) = (m*/M )a; c and at long wave lengths obeys 



c4(k) 



M 



2ne 2 
ea c M ' 



-k 



3^ + A. 
,— / 

A/n 



0(fc 3 ), 
(10) 



where 99 is the azimuthal angle of the wave vector k in the 
XY plane. The gapless mode on the other hand obeys 



e(m*w t 



-( M -Isin 2 2^)k 3 



+ 0(fc 4 ). (11) 



Within a low Zecman energy limit, where the mass of 
skyrmion is small, the gapped mode goes toward the 
higher energies and the effect of the mass of skyrmions 
becomes less significant. At this limit the prediction of 
the Chern-Simon theory, approaches to the current pre- 
diction of the microscopic Hartree-Fock approximation. 
However, for higher Zeeman energies, it is not clear if the 
prediction of the microscopic Hartree-Fock approxima- 
tion leads to a single gapless mode or it can support the 
gapped mode and subsequently the possibility of the ex- 
istence of a non-zero skyrmion mass too fl5|l . The square 
lattice can be unstable against the lattice fluctuations. 
Moreover the contribution of the exchange interaction to 
the Lame coefficients falls off exponentially (much faster 
than the contribution from the direct interaction) and 
when v — > 1 a triangular Wigner crystal configuration 
is more favorable than a square lattice. A similar situa- 
tion arises when either g — > or g — > 00. Our numerical 
studies show that this instability can be seen for values 
of the filling fraction and g-factors that are shown in Fig. 
||. The curve of R = 2A is also plotted. The overlap be- 
tween skyrmions is significant below this curve and the 
topological interaction becomes strong. The assumption 
of weakly interacting skyrmions may fail within the re- 
gion below this curve, and one should take the next to the 
linear topological XY terms into account. As it is shown 
in Fig. ||, for 1 — v\ < 0.01, and for certain direction 
of k, the gapless mode becomes imaginary, implying an 
instability of the square Skyrme lattice that is again con- 
sistent with the microscopic Hartree-Fock models [12|. 



The stable region of the square lattices is characterized 
by dark area and the stable region of the triangular lat- 



tices is shown by the white area. 20 It is seen the trian- 
gular lattice reappears when the Zeeman energy is small 
enough. Without Coulomb interactions a lattice config- 
uration is always unstable against the attractive interac- 
tion between skyrmions, i.e., a single skyrmion with topo- 
logical charge N is the global minima of the skyrmionic 
energy functional. 
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FIG. 2. The phase diagram for massive (and massless) 
skyrmions is shown. The stable region of the square skyrmion 
lattices are represented by the points (dark region). The tri- 
angular lattice is stable in the white region. For 1 — v\ < 0.01, 
the square Skyrme lattices are unstable. The curve of R — 2A 
is shown. The overlap between skyrmions is significant below 
this curve and the topological interaction becomes strong. 

As it is seen in Fig. [l], the quantum Hall skyrmions 
are massless at g = 0. For the massless case the gapped 
mode goes to infinity and there is just only one gapless 
mode 



(mu; c w(k)) z 



/27re 2 



V ek 

=2 



+ a 2 [A + 2/z] k 



-27 f + al\fi + A + -7] J k 2 x k 2 y + 0(fc 6 ), (12) 

which is identical to the gapless mode of the massive 
skyrmions. In fact, the long wave length power-law be- 
havior of the gapless mode is totally unaffected by the 
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mass of skyrmions. For both the massless and massive 
theories, we thus find that w cx fc 3 / 2 at long wave lengths, 
which is consistent with the microscopic Hartree-Fock 
calculations |l2| . Experimentally, the mass of skyrmions 
may not be probed directly by phonon excitations unless 
we excite the gapped modes (the cyclotron modes), e.g., 
by optical measurements JTsj ] . Furthermore, the mass of 
skyrmions can not change the melting point of the 2D 
skyrmions, since T m is specified by the elastic constants 
of Skyrme lattices (the 2D melting point just depends on 
the interaction energy between skyrmions) . Returning to 
the evaluation of the zero-point energy of the phonons, 
the dominant contribution at low temperatures is thus 
obtained as Er u = J^ fc ?iu;_(k). Similarly the average 
square displacement of the skyrmions is given by 



(u 2 ) 



-DxxO^) H~ Dyy 



(k) 



keBZ 



x [2n B (/?fiw-(k)) + l] 



(13) 



where ub{x) is the Bose-Einstein distribution function. 



B. Magnons 

Next we consider the spin waves of the Skyrme lattices 
by expanding the energy functional in Eq. (|J) also in the 
orientation field 6. Up to the quadratic order, it gives 
the XY-energy contribution 



E[u,9] = E[u] + \k xy J d 2 r|V0(r)| 2 , 



(14) 



where Kxy — J(kR) is the effective stiffness associated 
with gradients in the skyrmion orientations. As is shown 
in Ref. || , the effective action for the spin waves contain 
also a mass term (Ao is the moment of inertia) and we 
obtain finally 



s cS m 



f n /AM 

dr / d 2 r d T 

n J \ flc 



2a a 



K 



XY 



(15) 



The first term in Eq.([l5|) is the usual (dynamical) Berry's 
phase of a quantum Hall ferromagnet, where AM is the 
average change in the total magnetization induces by a 
single skyrmion texture. Note that in principle there is 
also a contribution from the Hopf term in the effective 
action for the quantum Hall ferromagnet, which at the 
quantum level ensures that the skyrmion obeys the cor- 
rect spin-statistics relation |l^]. Since both these terms 
are a total derivative, however, the equation of motion 
is not affected by these terms and the long wave length 
dispersion relation that follows from action ([l5|) turns 
out to be uj(k) = c s k, where c s = ^a c KxY/^o is the 



velocity of the spin waves. The contribution from these 
fluctuations to the total energy of the crystal is again 
EftQ = J2 k 7iw(k), and the mean square value of the as- 
sociated fluctuations is 



2 > = ^£ 



i 



keBZ 



A w(k) 



[2n s (/3Mk)) + 1] 



(16) 



The coupling between the displacement and the orien- 
tational fields (u and 9), turn out to be the next to 
leading order terms and are therefore negligible for our 
purposes. These terms lead to interactions between the 
phonons and spin waves, which we do not consider here. 
As a result we find for the zero-point energy of the 
phonons and spin waves simply Eji — Efi u + Efig. ^From 
Eq. ((l^), we find (u 2 ) ~ |1 — ^|i? 2 /5 at zero temperature. 
We also find (6 2 ) = ^/2U/K X y horn Eq. ^ where 
U = l/(2Ao). For small values of Kxy and/or large 
values of U the quantum fluctuations are severe that 
the disordered phase can emerge. For very small Zee- 
man energy, where a phase transition from square single- 
skyrmion lattice into a triangular bi-skyrmion lattice can 
be observed Jl|,§, we obtain (6 2 ) ~ g 0A6 . At the limit 
of small <?, fluctuations are negligible and the long-range 
order of the bi-skyrmion lattices is not influenced by the 
quantum fluctuations. 



IV. CONCLUSION 

In this paper we have studied a system of 2-dimensional 
quantum Hall skyrmions, starting from a Chern-Simon- 
Landau-Ginsburg mean field theory. A Maxwell term can 
be obtained through the gradient expansion of the Chern- 
Simon action around its minimum energy solution. This 
term is responsible for generating of the skyrmionic incr- 
tial mass. Away from v = 1 skyrmions stay in a crystal 
form. The long range order of crystal depends on the 
Landau level filling factor, and the Zeeman energy. Op- 
tical phonons which are out of phase fluctuations of the 
skyrmion lattices, are gapped since skyrmions carry an 
inertial mass. Therefore quantum Hall skyrmions be- 
have like a Wigner crystal in the presence of external 
magnetic field. The inertial mass of skyrmions is van- 
ished at zero Zeeman energy. In this situation the op- 
tical phonons are highly gapped, and they become in- 
accessible. We finish this paper with a final comment. 
At T = 0, and far from u = 1, the zero point quantum 
fluctuations of the phonons destroy the long-range or- 
der of skyrmion crystals. This is a high density limit of 
skyrmions where the Coulomb interaction screens from 
one 1/r to ln(l/r) by the Chern-Simon fluctuations. In 
this limit, the skyrmions behave like a gas of classical 
particles, i.e., they are crystalized under high pressures. 
This yields to the possibility of observing the re-entrance 
of the solid phase, followed by the disorder (liquid) phase 
at T = when v is far enough from 1. 
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